Enhanced ionization in small rare gas clusters 
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A detailed theoretical investigation of rare gas atom clusters under intense short laser pulses 
reveals that the mechanism of energy absorption is akin to enhanced ionization first discovered for 
diatomic molecules. The phenomenon is robust under changes of the atomic element (neon, argon, 
krypton, xenon), the number of atoms in the cluster (16 to 30 atoms have been studied) and the 
fluency of the laser pulse. In contrast to molecules it does not dissappear for circular polarization. 
We develop an analytical model relating the pulse length for maximum ionization to characteristic 
parameters of the cluster. 
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I. INTRODUCTION 

Building the bridge between atomic and solid state 
physics, cluster physics has become a vivid research field 
of its own. While the static properties of clusters are 
by now well understood, there remain many open prob- 
lems concerning the dynamics of clusters under external 
perturbations. In the case of weak perturbations, linear 
response theory has proven to be a valid tool for the in- 
vestigation of dynamical properties JjJ ; but with increas- 
ing strength of the perturbation the description of the 
cluster evolution becomes more and more involved || . 

On the other hand, experimental studies of (mostly) 
rare gas cluster interaction with highly charged projec- 
tiles as well as with short, intense laser pulses have 
produced a number of interesting results calling for an 
explanation. Experiments exploring the interaction of 
rare gas clusters with intense laser light have shown a 
big increase of energy absorption compared to the sin- 
gle atom case ||, |[ ||, ||, 0, g|. When irradiated with 
a 10 15 W/cm 2 femtosecond laser whose wavelength is in 
the optical regime, one observes, depending on the clus- 
ter size and the atomic element, ionic charge states of up 
to 40. These high charge states let the fragmenting ions 
gain an enormous amount of kinetic energy. The most 
spectacular example for this highly energetic process has 
certainly been the recent experimental observation of nu- 
clear fusion in a cluster ||. 

Here, we focus on clusters of some 10 atoms. We have 
developed a model containing the essential features of the 
interaction between the cluster and the laser field. It will 
be shown that energy absorption from the laser pulse pro- 
ceeds through a mechanism originally discovered for di- 
atomic molecules (enhanced ionization (ENIO) ]l0| , p"T|), 
whose generalization for the case of small clusters will be 
presented. 

The paper is organized as follows: after introducing 
the numerical model and comparing it to other types of 
simulations already existing in the literature in section II, 
we investigate the dependence of energy absorption and 
ionization yield on the pulse length in a series of clusters 
in section III. From the results of these calculations a 
generic behavior emerges which can be explained by in- 



voking the above-mentioned enhanced ionization mecha- 
nism as explained in section IV. We give strong evidence 
that this mechanism should play an important role in the 
laser-cluster interaction over a large range of parameters 
detailed in section V. Finally we condense our picture of 
the ionization process into a simple analytical expression 
which quantifies the role of the experimentally accessible 
variables like cluster size or atomic element on the pro- 
cess of energy absorption in section VI. The last section 
VII summarizes our work. Atomic units are used if not 
stated otherwise. 



II. THE CLUSTER MODEL 

A. Theoretical formulation and numerical 
implementation 

Since the dimension of the problem is far too high to al- 
low for an exact quantum mechanical treatment we have 
formulated a model to describe the dynamics of rare gas 
clusters in strong laser fields. We resort to a classical 
treatment, with a few but essential quantum mechanical 
elements. 

Initially, before the onset of the laser pulse, we find 
the equilibrium ground state configuration of the clus- 
ter with Lennard- Jones interaction between the neutral 
cluster atoms. The global minima for this potential are 
readily available. The electrons are assumed to be local- 
ized at the nuclei. 

After fixing the initial shape of the cluster, we start 
the time evolution switching on the laser pulse |l^ ]. For 
the electrons, this evolution consists of two parts: firstly, 
the modeling of the bound state and the process of ion- 
ization from this state; secondly, the propagation after 
being ionized from an atom. We will refer to the first 
process as inner ionization, in contrast to the outer ion- 
ization, which has the effect that an electron leaves the 
cluster. Inner ionization contains processes beyond clas- 
sical mechanics, while the subsequent propagation and 
(eventually) outer ionization is described classically via 
integration of Newton's equations. 

When irradiating a cluster with intense laser light, two 
processes can lead, at least in principle, to inner ioniza- 
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tion: field ionization and electron impact ionization. In 
the case of field ionization, the electric field inside the 
cluster (initially only the laser field, later the combined 
field of laser, ions and electrons) leads to a lowering of 
the potential barriers, so that an electron can leave its 
mother atom via tunneling |fl3| or even over-the-barrier 
Jj"4| |. Electrons which are already inner ionized, but not 
yet outer ionized, can further lead to electron impact 
ionization. This mechanism was shown to play almost 
no role in small clusters |l5|, as the average free path 
length with respect to electron impact ionization is much 
larger than the cluster radius. For this reason we only 
consider field ionization in our model. 

The model is implemented as follows. Before the onset 
of the pulse, the electrons of the cluster are assumed to be 
localized at the atomic positions; the bare nucleus and all 
the electrons of an atom are treated as one neutral clas- 
sical particle. It is only later that the electrons (through 
inner ionization) are born as separate classical particles. 
Hence, the number of particles in our simulation changes 
with time. 

Two classical charged particles at positions r\ and f% 
interact via the potential 



V{r\,rz) = Z\Z 2 (\fi - r 2 \ 2 + a Zl + a z ^j 



-1/2 



(1) 



where Z\ and Z 2 are the charges of the two particles. The 
azi arc softcore parameters, which help to regularize the 
Coulomb singularity. We are using a Z-dependent a, so 
that the depth of the atomic potential can be adjusted to 
the current binding energy. To determine a^, we proceed 
as follows: 

- For electrons, a has been chosen to be 0.1, i.e. 
o_i = 0.1 

- For an electron which is localized in the minimum of 
the potential of an atom with charge Z, we demand 
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with a small positive parameter e = 0.01. Then we 
can solve for az 



az 



Kind + ^bind + e 2 



a-i 



(3) 



At every time step dt, we calculate for each atom the 
ionization probability of the outermost electron in the 
following fashion: if Bj is the total electric field at the 
position Rj of atom j, 
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FIG. 1: Tunneling from a bound state with energy E n in the 
potential V(r) 



with the polarization vector e p of the laser field, we cal- 
culate the tunneling integral in the direction of B: 



I(t) = exp 
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y/2(V(r) - E n ) dr 



(5) 



The role of n and r 2 is shown in Fig. [l] where V(r) con- 
tains the Coulombic potential terms and those coming 
from the laser field. The energy level E n is defined as 
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Hence, the atomic energy level is shifted by the surround- 
ing charges and the laser potential; the potential of the 
atom from which the electron will be ionized has to be 
subtracted because its influence is already contained in 
E^ tora . The positions r\ and r 2 , where E n crosses the 
potential curve V(r), are determined numerically; the 
search for r 2 is continued until I(t) < 10~ 10 . If we find a 
position r' along the direction of B with dV/dr\ r=r > = 
and E n > V(r) for r < r', over-the-barrier ionization is 
possible. In this case we simply put I(t) = 1. The tunnel- 
ing rate w(t) is then the tunneling probability I(t) mul- 
tiplied with the frequency of the electron hitting the po- 
tential barrier. In a semiclassical picture, this frequency 
is just the inverse of the Kepler period T„ belonging to 
an orbit with binding energy E n : 



T n = 7r(ZAtom 4 
Hence, the tunneling rate is 
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w(t) = I(t)/T n . 



(7) 



(8) 



The tunneling probability over a unit of time dt is 
P{t) = w(t) dt. By comparison with a random number 
z (is P{t) > zl) we decide if the electron in question 
tunnels. If so, we place the electron, which now becomes 
a classical particle, outside the potential barrier as close 
to r 2 as possible, with the exact position and momentum 
of the electron determined by conservation of the total 
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FIG. 2: Number of classically treated particles (a) and total 
charge of the cluster (b). 



energy |16| . If the ionization happens to be over-the- 
barrier, we put the electron on top of the barrier, where 
dV/dr = 0. The atomic charge is raised by 1 and the 
next virtual electron is allowed to tunnel. 

The particles are classically propagated by integrating 
Newton's equations of motion. We have used a symplec- 
tic integrator |l7| with a time step of dt = 0.1. 

B. A typical run 

Although later we will use a Monte Carlo ensemble 
to calculate experimentally accessible observables, for a 
qualitative understanding of the phenomena it is suffi- 
cient to have a closer look on a single event, since the 
overall behavior of the ensemble members is quite sim- 
ilar. As an example we consider a Nei6-cluster. The 
applied pulse has a peak intensity of I = 10 15 W/cm 2 , a 
frequency of u> = 0.055 a.u. (780 nm) and it extends over 
20 cycles, so that the pulse length is approx. T w 55 fs. 
We chose a sin 2 -function for the pulse envelope, i.e. the 
pulse is of the form 

/(f) = F sin 2 (nt/T) sin(wf) (0 < f < T). 

(9) 

Fig. H shows the total number of classically treated par- 
ticles and the cluster charge, i.e. the number of electrons 



which have left the cluster, as a function of time. After 
approximately 750 fs the intensity of the laser is suffi- 
ciently high for the first inner ionization event, followed 
by a rapid increase of the number of classical particles as 
well as the cluster charge. Obviously, the ionization of 
the first few electrons leads to an 'avalanche effect': the 
inner ionized electrons create a strong electric field inside 
the cluster, which helps to inner ionize further electrons 
(this is reminiscent of the ionization ignition mechanism 

H). 

Fig. U shows the energy absorbed by the cluster and 
the mean interionic distance, again as a function of time. 
If Jtf is the set of nuclei with mass M, S the set of inner 
ionized electrons (mass m) already treated classically and 

the set of electrons which are still bound, the cluster 
energy is defined as 

y^-+y^- + y vu+ye?** no) 

^ 2M ^ 2m ^ lJ ^ 1 v ; 

+ J2 z i (?p ■ Ri) fit) - E & ■ fit) • (n) 

The absorbed energy is the difference of the total energies 
before and after the laser pulse. As we can see, this rather 
small cluster can already absorb a considerable amount 
of energy. The oscillations are due to the ponderomotive 
potential and have no direct influence on the net energy 
absorption. As the cluster gets charged, it begins to ex- 
pand, i.e. the mean interionic distance will increase. For 
a cluster consisting of N atoms, it is defined as 

/ N XV2 

^) = hv£?^ {l ^-^' |2} ) • (12) 

At the intensity used here, the cluster disintegrates 
completely, i.e. we observe only atomic fragments after 
the pulse. Note that the expansion of the cluster takes 
place adiabatically compared to the time scales of the 
laser frequency and the electronic motion, but on the 
same time scale as the pulse length. Hence, it is possible 
to explore radius-dependent properties of the cluster by 
varying the pulse length. 

C. Comparison to other models 

Similar models for the theoretical description of small 
rare gas clusters in strong laser fields have been proposed 
by other authors. The common feature is the use of clas- 
sical mechanics as the basic ingredient of the description; 
the main differences lie in the treatment of inner ioniza- 
tion. 



1. The onion model 

Following its name this model fl8|| , as ours, treats the 
weakest bound electron per atom or ion as the active 
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FIG. 3: Absorbed energy (a) and mean interionic distance (b) 
(see Eq. (Jl2|)) as a function of time. 



electron which can be inner ionized. However, tunnel- 
ing is not included. Instead, the active electron circles 
its mother ion on a Kepler orbit with an energy cor- 
responding to the quantum mechanical binding energy. 
This orbit is then deformed in time by the laser field and 
the surrounding charges, and inner ionization can take 
place, which is defined by reaching a certain distance r c 
from the mother ion. If the electron exceeds this dis- 
tance, it is assumed to be inner ionized, and the next 
electron is put onto a (now deeper bound) Kepler orbit. 

In common with our model, inner ionization takes 
place sequentially, and the already inner ionized electrons 
are propagated classically. By neglecting any tunneling 
contributions, the first inner ionization will take place a 
certain time St later than in our case. As we will see, 
this delay can have a significant influence on the subse- 
quent dynamics of the cluster rendering a fully classical 
treatment as in the onion model problematic. 



2. Tunneling via the Landau rate 

In this ansatz p9[ the tunneling mechanism is taken 
into account, but the tunneling integral is not calculated 
explicitely. Instead, the time dependent tunneling prob- 



ability is estimated using the Landau tunneling rate p0| , 
which has been established for the case of an atom in a 
quasistatic laser field. For this purpose, the total electric 
field (laser + electrons + ions) at the position of an atom 
is used. For inner ionization to happen in this model the 
field strength has to be strong enough at just one point, 
namely, the position of the atom. Hence, an electron 
which comes by chance close to an ion will create locally 
such a strong field that ionization can hardly be avoided. 
In our model, in contrast, the entire environment of an 
atom must be suitable for ionization. Consequently, we 
get lower ionization rates than in |l9[] but in better agree- 
ment with Jl5| . 

3. Using the full Coulomb potential 

In a third model |H| , the full tunneling integral is cal- 
culated at each time step as in our case. Moreover, the 
authors do not use a softcore potential but employ the 
full Coulomb potential and regularize the singular equa- 
tions of motion. The use of the full Coulomb potential 
may, at first glance, seem to be an advantage. However, 
it leads to the problem of classical autoionization: a clas- 
sical electron in a Coulomb potential can be arbitrarily 
deeply bound, since there is no uncertainty relation. The 
corresponding energy gained is available as kinetic energy 
for the ionization of other electrons. In our model, this 
unphysical behavior is avoided by requiring the minimum 
of the softcore potential to coincide with the quantum 
mechanical binding energy. 

In fl5|| a recapture mechanism was built in to account 
for this problem, which can transform already classically 
treated electrons back into a virtual bound state exis- 
tence. Not to mention the peculiarities which arise when 
defining the exact conditions for recapture in a many- 
particle environment, the numerical expenditure is higher 
than in our case. Nevertheless, the ionization yields cal- 
culated with this model agree quite well with our results. 

III. ABSORPTION PROPERTIES FOR 
DIFFERENT PULSE LENGTHS 

To investigate how the expansion of a cluster during 
the interaction with a strong laser pulse influences its ab- 
sorption behavior, we have calculated the absorbed en- 
ergy and the average ionic charges after the interaction 
for various pulse lengths. To keep the amount of energy 
delivered by these pulses fixed, we demand the fluency to 
be constant, i.e. 

E(T):= [ f(t) dt = const. (13) 
Jo 

For a pulse of the shape of Eq. (||), we obtain 

E(t) = 3F 2 T/16. (14) 
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FIG. 4: Energy absorption of Nei6 (°), Ari6 (•), Kris (o), 
and Xei6 (*) for different pulse lengths (see text). 




FIG. 5: Average atom charge as a function of pulse lengths 
for different clusters, see Fig. ^. The solid lines are to guide 
the eye. 



atoms. This changes for heavier rare gases. 



B. Generic cluster behavior in Arie 

The absorbed energy (Fig. ||) and average atomic 
charge (Fig. ||) for an Ari6 cluster decrease monotoni- 
cally for very short pulse lengths only. Both quantities 
start to rise again for T w 2000 a.u., reaching a maxi- 
mum at an optimal pulse length of T sa 4500 a.u. (we will 
refer to this optimal pulse length as T* subsequently) be- 
fore they ultimately decrease. The average charge state 
reached at T = T* is around 4.5, which is considerably 
higher than what one would expect from a single Argon 
atom. 

The maximum is much more pronounced for the ion- 
ization yield than for the absorbed energy for two rea- 
sons. Firstly, for shorter pulses the mean internuclear 
distance just after the pulse will in general be smaller 
than for longer pulses, so that the Coulomb explosion 
energy will increase with decreasing pulse length. Sec- 
ondly, the ionized electrons, which can be considered to 
be quasi-free, acquire a higher kinetic energy for shorter 
pulses, since the intensity is higher than for longer pulses. 
These two effects wash out the minimum in the curve for 
the absorbed energy and decrease the contrast between 
minimum and maximum. 

We will postpone the further discussion and explana- 
tion of this structure after we have taken a look at the 
corresponding results for Krig and Xei6- 



Same qualitative behavior as for Argon: Kris 
and Xei6 



For the reference pulse, we chose the parameters already 
used in the single run from the previous section: F = 0.16 
a.u., u> = 0.055 a.u. and a pulse length of 20 cycles. 
Shorter pulses have a higher maximum field strength and 
longer pulses a lower maximum field strength, respec- 
tively. The results were obtained by averaging over a 
Monte Carlo ensemble consisting of 20 clusters. 



A. The light Nei6 cluster: almost atomic behavior 

The first cluster we consider is Nei6- The absorbed 
energy and the average atomic charges as a function of 
pulse length under the constraint of Eq. ([l3|) are shown 
in Figs. |] and |[ respectively. 

Both curves decrease monotonically with increasing 
pulse length. The average atomic charge exhibits a 
plateau for pulse lengths around T — 3000 a.u.. At this 
point, we can already suspect that the plateau might have 
its origin in the expansion of the cluster during the pulse. 
However, in the case of Nei6 its effect seems to be rather 
small. In fact, Neon clusters of this size behave more like 



Qualitatively, the behavior for Krig as well as for Xei6 
is the same as for Ari6: the energy absorption and the 
average atomic charge peak for a certain pulse length T* , 
which depends on the type of the cluster. The average 
atomic charge increases with the atomic mass, reaching 
an average of almost 7 + for the Xei6 cluster; this is of 
course due to the decreasing binding energies of the elec- 
trons in the heavier elements. The value of T* is chang- 
ing for different clusters; we will come back to this point 
later, after we have understood the reason for the exis- 
tence of T*. 



IV. CALCULATIONS WITH FIXED NUCLEI 

The existence of an optimal pulse length T* can be re- 
lated to an optimal cluster geometry with critical radius 
R* , which maximizes the energy absorption (and also the 
ionization) of a small rare gas cluster in a strong field. 
Given the existence of such a critical radius, which should 
be larger than the equilibrium radius Rq, the occurrence 
of T* can be readily explained. 
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FIG. 6: Cluster charge, calculated with fixed nuclei, as a 
function of the mean interionic distance (see Eq. (Jl^) , in units 
of the equilibrium mean interionic distance Ro) for Nei6 (o), 
Ane (•), Kn 6 (o) and Xei 6 (*). 



For very short pulse lengths, the cluster has almost 
no time to expand during the pulse, so that the critical 
radius will be reached only well after the pulse is already 
switched off. For longer and longer pulses, the cluster 
will reach R* at earlier and earlier times. For a certain 
pulse length T*, the time of reaching R* will roughly 
coincide with the maximum of the pulse, which leads 
to optimal absorption. If the pulses are becoming even 
longer, R* will be reached already before the maximum of 
the pulse and the maximum intensity is decreasing due to 
the energy normalization (Eq. (13)). Both effects lead to 
a decrease in energy absorption as well as in the average 
atomic ionization for T > T*. 

What remains to be shown is that the critical radius 
R* really exists and to explain its origin. To this end, we 
have calculated the cluster ionization yield for different, 
but fixed cluster radii. We have accomplished this by 
applying a scaling transformation 



R4 =>■ XRa 



(15) 



to the atomic positions, with A = 1 corresponding to the 
ground state configuration. The pulse we use is the ref- 
erence pulse for the calculations of the last section, i.e. 
of the form Eq. (|) with F = 0.16 a.u., u = 0.055 a.u. 
and 20 cycles length. The results of these calculations 
are shown in Fig. ^. For all four clusters under con- 
sideration, we observe the existence of a critical radius 
R* ^J, which is larger than the equilibrium one. This 
means that it will be reached during the expansion the 
clusters undergo when irradiated with intense laser light. 
The results of the previous section thus can indeed be ex- 
plained by the existence of R* . The question remaining 
to be settled is: which mechanism is responsible for R*l 
In principle, two possibilities are available: first, it could 
be a resonance effect where the electrons inside the clus- 
ter oscillate at a certain characteristic frequency which 
would coincide with the laser frequency at a certain clus- 
ter size. This kind of mechanism is well known from the 



plasmon resonance in metal clusters p2| ; being originally 
a weak-field-concept, the plasmon has been claimed to 
play an important role also in the strong field regime M . 
However, one needs a delocalized electron cloud to create 
a plasmon resonance; in rare gas clusters, this condition 
is not fulfilled. 

The second mechanism would be a generalization of 
a concept first discovered for linear, diatomic molecules 
]To| , pd] , [23] ] , called enhanced ionization (ENIO ) . It can be 
qualitatively explained by looking at the potential curve 
of such a molecule exposed to a quasistatic electric field 
(see Fig. Q): the upper energy level of the two levels 
1(7+ and lcr_, which emanate from the bonding and the 
antibonding molecular orbital when an electric field is 
switched on, will lie above the inner potential barrier but 
below the outer potential barrier when the internuclear 
distance R is rather small. On the other hand, when R 
is rather large, the level will lie below the inner barrier 
but above the outer barrier (using the terms introduced 
in the previous section, we can say that inner ionization 
is easier than outer ionization for small R and vice versa 
for large R). For an intermediate value of R, typically 
around 6-8 a.u., the interplay between inner and outer 
ionization will lead to a maximum in the ionization rate. 

This mechanism has been shown to be operative not 
only in linear molecules, but also in triatomic molecules 
of triangular shape [^4|, In this case, the simple pic- 
ture of Fig. |?] is already slightly distorted, and it is more 
appropriate to think of enhanced ionization in terms of 
an optimal balance between inner and outer ionization, 
which makes the generalization of the mechanism to a 
true many-body system like a cluster much easier. 

One characteristic feature of the enhanced ionization 
mechanism is its relative insensitivity on the frequency 
of the applied laser field. As long as the quasistatic 
picture is valid, the value of R* should not change sig- 
nificantly with the laser frequency. On the other hand, 
any resonance-type mechanism like the plasmon picture 
should exhibit a strong dependence of R* on the laser fre- 
quency. We have calculated the ionization yield of Nei6 
and Ari6 for three different frequencies, with the results 
shown in Fig. |[ The position of R* does hardly change 
for different laser frequencies. Hence, we can exclude any 
kind of resonance behavior in favor of the enhanced ion- 
ization mechanism. 

Although the position of R* does not change with the 
laser frequency, the ionization yield does. This is due 
to the fact that electrons which are already outer ionized 
tend to leave the cluster region faster when the frequency 
is smaller: the quiver amplitude of an electron in an elec- 
tric field of frequency u> is proportional to 1/cj 2 . Hence, 
on average in fields of higher frequencies the already ion- 
ized electrons will stay closer to the cluster for a longer 
time and lead to an increased field ionization rate. 
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FIG. 7: Schematic potential curves and the upper energy level 
1<t + of a diatomic molecule for different internuclear distances. 
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FIG. 8: Atomic ionization yield for the three frequencies u) = 
0.055 a.u. (solid line), ui = 0.075 a.u. (long dashed line) and 
Lu — 0.11 a.u. (dashed line) for Nei6 (a) and Ari6 (b). The 
pulse length was T = 55 fs. 



V. EXPLORATION OF THE PARAMETERS 
CONTROLLING LASER-CLUSTER 
INTERACTION 

Having established the basic mechanism for coupling 
energy from the laser pulse into small rare gas clusters, 



FIG. 9: Average atomic charge state as a function of the mean 
interionic distance for Nei6 (•), Ne2o (o), Ne25 (o) and Ne3o 
(*)• 



> 300 




FIG. 10: Absorbed energy/atom as a function of the mean 
interionic distance for Nei6 (•), Ne2o (°), Ne25 (o) and Ne3o 
(*)■ 



we will explore now the influence of different parameters 
on this mechanism, such as cluster size, energy content 
of the laser pulse and laser polarization. 



A. Different cluster sizes 

As in the previous section, we present the cluster re- 
sponse with fixed nuclei first and relate the results of 
these calculations to the absorption behavior when the 
nuclei are allowed to move. The average atomic charge 
and the absorbed energy were calculated as a function of 
the mean interionic distance again, analogous to section 
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The equilibrium value Rq does hardly change when 
going to bigger clusters. The variation of Ro for Nei6, 
Ne2o, Ne25 an d Ne3o is only about 0.01 a.u.. We have 
used the same pulse as in section IV. As can be seen 
from Fig. |9| and Fig. the bigger clusters show almost 
no difference compared to Nei6 when the observables are 
normalized on the number of cluster atoms. In particular 
the existence of a critical distance R* > Ro is confirmed 
in all cases. 

There is no hint on a transition to a collective behavior 
at these cluster sizes. If we think of cluster physics as the 
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transition regime between atomic and solid state physics, 
we are still on the atomic side with a cluster of 30 atoms. 

From the fact that the charge per atom is almost in- 
dependent of the number of cluster atoms we may con- 
clude that only the next-neighbor-atoms participate in 
the mechanism of enhanced ionization; otherwise the ef- 
fectivity of this mechanism should change with the clus- 
ter size. The absorbed energy per atom, however, is vary- 
ing with the number of atoms. This effect can be easily 
explained by calculating the change in the potential en- 
ergy U (N) of a cluster consisting of N ions of charge Z 
and radius R if one adds a new ion with the same charge 
Z to the cluster. If one assumes that this new ion is 
placed at the border of the cluster, then 



U(N + 1) = U(N) + NZ 2 /R. 



(16) 



If 47rri? denotes the volume per atom, then R — N 1 ^! 
and 



U{N + 1) = U(N) + N 2/3 Z 2 /r s 



(17) 



With A as a continuous variable one is left with the dif- 
ferential equation 



dU(N) _ N 2 /' 3 Z 2 



so that finally 



dN 



U{N) . IH^E. 

5 r s 



(18) 



(19) 



Hence, the potential energy per atom U/N has increases 
with N 2 / z if the charge per atom is independent of TV. 

Proceeding from Ne to Ar clusters, one finds again 
that the effectivity of the ionization mechanism hardly 
changes when changing the cluster size, while the ab- 
sorbed energy per atom increases with N, for the same 
reason as discussed above (see Fig. [ll] and Fig. [l2|) . How- 
ever, while the Ne clusters show only a little shift of R* as 
a function of cluster size, the ratio of R* to Rq increases 
slightly more with increasing N for Ar. This is probably 
due to a larger downshift of the atomic energy levels by 
the increased total amount of surrounding charge when 
N is increased. As can be seen from Fig. [?], a downshift 
of the atomic energy levels leads to an increase in R*. 
Since the electron release in argon clusters is larger than 
in Ne clusters higher charged ions are generated than in 
neon clusters rendering this effect more pronounced for 
Ar clusters. 

Since we have found a critical cluster radius R* with 
R* > Ro in all cases considered, it is not too surprising 
that we find a behavior analogous to the small clusters 
of Fig. |] and Fig. || if the bigger clusters are allowed 
to expand freely. The results of these calculations, with 
the pulse normalization being identical to section |ll|, are 
shown in Fig. |l^ and Fig. [l4|. For clarity, we have plotted 
the total cluster charge instead of the average atomic 
charge, which is almost the same independent of N. 




FIG. 11: Average atomic charge state as a function of the 
mean interionic distance for Arie (•), Ar2o (°), Ar25 (o) and 
Ar 30 (*). 




R/Ro 



FIG. 12: Absorbed energy /atom as a function of the mean 
interionic distance for Ari6 (•), Ar2o (°), Ar25 (o) and Ar3o 
(*)• 



The overall structure of the curves is seen to be quite 
similar throughout the different cluster sizes, as one 
would expect from Fig. | and Fig. [ll]. In the case of Ne 
clusters, the plateau which has been observed for Nei6 
already in Fig. ^ goes over into a small maximum with 
increasing A, which indicates that the enhanced ioniza- 
tion mechanism is slightly more efficient for larger clus- 
ters when the ions are allowed to move. One tendency 
which can be observed for the Ar clusters is that T* in- 
creases with increasing N. We have seen in Fig. [ll] that 
R* increases also with N, so that the larger Ar clusters 
have to travel a longer distance until they reach the crit- 
ical radius. For the Ne clusters, on the contrary, the 
curves show almost no shift in the T-direction when N 
is changed. We will investigate the dependence of the 
expansion process on the various cluster parameters like 



size and atom charge in closer detail in section VI 



B. The influence of the pulse normalization 



Changing the laser intensity / in the case of H^" 



just a single electron available, leads to 
R* when / is increased and vice versa fl2l 



with 
a decrease of 
. For clusters 
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FIG. 14: Cluster charge as a function of pulse length for Ari6 
(•), Ar 2 o (°), Ar 2 s (o) and Ar 30 (*). Lines are to guide the 
eye. 



the situation is much more complicated because with in- 
creasing / lower lying energy levels will be ionized, so 
that it is a priori not clear in which way a change of 
the laser intensity (in a calculation with fixed nuclei) will 
influence the value of R*. Fig. [l5| shows the static ion- 
ization yields for Neie, Ari6, Kri6 and Xei6 under the 
influence of the pulse used so far (i.e. a peak intensity 
of I\ = 8.99 • 10 14 W/cm ), compared to the result of a 
calculation with I 2 = 2.19- 10 15 W/cm (in both cases the 
pulse was of the form Eq. (||) with lo = 0.055 a.u. and 
T = 55 fs). In all four cases R* is larger than Rq and 
can be reached by cluster expansion. The value of R* 
is, if at all, only slightly decreased in the case of higher 
intensity: due to the large number of electrons involved 
the geometry of the problem is obviously not as sensitive 
to the laser field strength as in the Hj" case. 

Of course, the ionization yield is higher when the inten- 
sity is increased. This leads to significantly shorter ex- 
pansion times when the nuclei are allowed to move. Con- 
sequently, the optimal pulse lengths T* are now shifted 
towards smaller values, as can be seen in Fig. Hq in ac- 
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FIG. 15: Static ionization yield at the two intensities I\ = 
8.99 ■ 10 14 W/cm 2 («)andI 2 = 2.19 ■ 10 15 W/cm 2 (o) for Nei 6 
(a), Ari6 (b), Kri6 (c) and Xei6 (d). Lines are to guide the 
eye. 
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FIG. 16: Pulse length dependent ionization yields with a pulse 
energy corresponding to I2 = 2.19- 10 15 and a pulse length of 
20 cycles for Nei6 (o), Ari6 (•), Kri6 (o) and Xei6 (*). 



cordance with our picture of the ionization process. 



C. Enhanced ionization and circular polarization 

So far all the results presented are expected to hold 
also for diatomic molecules. One main difference between 
such a molecule and a cluster is the molecular axis: the 
whole picture of ENIO as sketched in Fig. [?| relies on the 
fact that the polarization direction of the applied laser 
field coincides with the internuclear axis. And indeed, 
experiments as well as calculations with a polarization 
axis perpendicular to the molecular axis have shown no 
signature of enhanced ionization |24|, ^?J. For the same 
reason ENIO is much less efficient under circular polar- 
ization. 

On the other hand, a cluster is (in first approximation) 
spherically symmetric. Thus one would expect enhanced 
ionization to work also with circularly polarized light. To 
test this hypothesis, we have performed the same calcula- 
tions as in the previous sections, but now with circularly 
polarized laser light. We have chosen the field strength 
of the laser such that the energy content of a pulse with 
a certain pulse length T remains constant when passing 
from linear to circular polarization: if 



F sin 2 '{ir/Tt) sin(o;i)e a 



(20) 



is a laser pulse with linear polarization in ^-direction, 
then 

F /V2 sin 2 (w/Tt) (sin(ut)e x + co$(wt)e v ) (21) 

is the corresponding circularly polarized pulse. With this 
definition the maximum field strength is decreased by a 
factor of As expected, ENIO also exists for circu- 
larly polarized laser pulses. Fig. |l7] shows the calcula- 
tions with fixed nuclei, Fig. [IJ the corresponding results 
with moving nuclei. In the case of static nuclei, we find 
that the ionization yield in the critical regime is almost 
as high for circular as for linear polarization (Fig. |^), 
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FIG. 17: Atomic charges with fixed nuclei and circular po- 
larization for Nei6 (o), Ari6 (•), Kri6 (o) and Xeie (*). The 
pulse parameters are Fo = 0.16/v2 a.u., U) = 0.055 a.u. and 
T = 55 fs. 
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FIG. 18: Same as Fig. ^ but with moving nuclei. 



in sharp contrast to the above-mentioned results for di- 
atomic molecules. Consequently, when the nuclei are 
allowed to move we also get qualitatively the same re- 
sults (Fig. |l8| ) as with linear polarization (Fig. [5]). It 
is only for rather long pulses that in the case of Nei6 
and Ari6 the ionization yield is significantly lower than 
in the linear case, which is due to the reduced maximum 
field strength. Summarizing our exploration of different 
parameters we find that ENIO for clusters is a rather ro- 
bust phenomenon. This has motivated us to ask if the 
optimum pulse length T* can be quantitatively linked to 
the critical radius R*. 



VI. ANALYTICAL FORMULA FOR THE 
COULOMB EXPLOSION 

To isolate the relation of T* to R* we divide the time- 
dependent dynamics into three different phases: phase I 
denotes the time from the onset of the laser pulse until 
50% of the atoms in the cluster have lost one electron due 
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FIG. 19: Sketch of phases I, II and III during the pulse (see 
text) 



to inner ionization. We will refer to this time as To sub- 
sequently. Since some of the inner ionized electrons will 
leave the cluster, we can say that To marks the beginning 
of the expansion process. 

In this first phase inner ionization is dominated by 
atomic processes, the environment plays only a minor 
role. For a single atom/ion the time-dependent proba- 
bility that the active electron is not yet ionized reads, in 
terms of the field- and binding energy dependent ioniza- 
tion rate w(f(t),Eb) 



neutral 



(t) — exp 



w(f(t'),E b ) dt' 



(22) 



where Ef, denotes the binding energy. The probability 
that in a cluster consisting of N such atoms no electron 
has been ionized is given by 



^neutral (*) ~ [-PneutralW] 



A? 



(23) 



The exponential dependence on N renders ^eutrSC*) 
practically a step function. Hence, the exact value (be- 
tween zero and one) for the definition of To is not rele- 
vant. We determine T from P^^(T ) = 1/2 which is 
tantamount to demanding that on average 50 % of the 
atoms in the cluster are singly ionized at To. 

The second phase contains the cluster expansion up to 
the critical time T*, when the critical cluster radius R* 
is reached. Hence, the critical time is the sum of To and 
the expansion time T exp : 



T* = Tn 



L cxp ■ 



(24) 



The third phase is finally the time from reaching R* to 
the end of the pulse. A schematic picture of the different 
phases is shown in Fig. 

From T) until the end of the pulse, the total cluster 
charge increases from Z = N/2 to Z = Z& na \. As a 
first approximation, we assume that the expansion from 
R = i?o to R = R* is driven by an effective charge 
per atom Z — aZg na i/JV with a constant factor a which 
stands for the efficiency of the enhanced ionization mech- 
anism. Furthermore, the expansion is assumed to be ac- 
complished by the Coulomb repulsion of the nuclei only, 



i.e. we neglect the influence of the laser field as well as of 
the electronic dynamics on the expansion process. Under 
these two assumptions, we can use energy conservation 
to write 



N 

E 

i=l 



M 



N 



z 2 



^ R 

(i#)=l 



«(*) 



E, 



(25) 



where Rij{t) = \Ri{t) — Rj(t)\, M is the atomic mass 
and Vi the respective atomic velocities. As a further ap- 
proximation we assume that the expansion takes place 
in a homogenous and isotropic way, so that it can be 
described by a common expansion parameter X(t) with 
Ri(t) = X(t)R t (0). Defining 

N 



K :=^21mR 2 (0) 



i=l 
N 



0¥i)=i 



(^final/A) 2 

%(0) 



(26) 



and taking into account that the kinetic energy is zero 
before the expansion we may write the energy balance of 
Eq. (H) as 



K X 2 (t) 



a 2 a 2 

W) Vo = W) Va 



(27) 



Finally, Eq. ( |27| ) may be rearranged as a differential equa- 
tion for A(i) 



d\{t) 
dt 



X(t)J K Q 



1/2 



(28) 



which can be solved analytically by separation of vari- 
ables for the expansion time: 



T- 




log (Va 



1 



x = l 



(29) 



The ratio Kq/Vq determines the time scale for the ex- 
pansion of the cluster. By replacing Ri 3 -(0), the distance 
between two ions in Vq, with the cluster radius R (which 
would be an exact approximation if all ions were placed 
on the surface of the cluster), we can estimate how this 
time scale depends on the characteristic variables of a 
cluster: 



MR A 



V {N-l){Z^/Nf 



(30) 



From this equation we can read off how the expan- 
sion process changes when the number of atoms N is 
changed while keeping all other parameters fixed: if 
Vatom = 4/37rrg is the volume of one atom inside the 
cluster, then R 3 — Nr 3 Hence, the time scale of the 
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FIG. 20: Expansion time (numerical data) as a function of 
(Ko/Vo) 1 ^ 2 /(A) and linear fits (see text). Two different en- 
ergy normalizations were used: Fq = 0.16 a.u. (solid line) 
and Fo — 0.25 a.u. (dashed line), both at a frequency of 
uj — 0.055 a.u. and a pulse length of T = 55 fs. Aris(«), 
Ar 20 (o), Ar 25 (D), Ar 30 (o), Kn 6 (*) and Xei 6 (x) 



pulse. This result a posteriori justifies the approxima- 
tions we have made in establishing our expansion model. 

The fact that a remains almost the same when chang- 
ing the pulse normalization is certainly an unexpected 
result; it is probably valid only for a limited range of 
pulse energy contents, if one thinks of a as a measure 
for the efficiency of ENIO. At least in the limit of a very 
large pulse energy, when the electric field of the laser 
is larger than the electric field from the charges in the 
cluster, we expect the ENIO mechanism to play no im- 
portant role any more, since the cluster geometry will be 
washed out. However, the good agreement of the linear 
fit in Fig. |o| with our numerical data for each of the two 
normalizations separately points to a deeper scaling re- 
lation between the various clusters, the reason for which 
will be explored in future work. 



VII. SUMMARY AND CONCLUSION 



expansion is governed by the factor N/(N — 1), which 
depends only weakly on N. This leads us to the conclu- 
sion that indeed the increase of T* with A" in Fig. |l4| is 
not due to the expansion process but due to the increase 
of R* with N only (see Fig. [□]). 

A. Scaling of the optimal pulse length 

With Eq. (H^) at hand we are able to set up a relation 
between the optimal pulse lengths for various cluster if 
we make one last assumption: the factor a, which de- 
termines the ratio between the average atomic charge Z 
driving the expansion up to R* and the final charge per 
atom after the pulse, Zs na \/N, is identical for all clusters. 
If this hypothesis was true, then 1 — a would be a univer- 
sal measure for the efficiency of the ENIO mechanism. 

If a is the same for all clusters, we get from Eq. ( |29| ) 
a linear relation between the expansion time T cxp and 
/(A)(To/Vb) 1 / 2 , different for each cluster. This predic- 
tion is confirmed by Fig. ^ which shows the expansion 
times T exp = T*/2 — To as a function of the cluster- 
dependent values of {Kq/Vq) 1 / 2 /(A) for different clusters. 
We have obtained A from the respective static calcula- 
tions for each cluster. A linear fit to the data yields 
a = 0.38 and a = 0.37 for energy normalized pulses at 
Fo = 0.16 and Fo = 0.25, respectively. The correlation 
coefficient is in both cases higher than 0.99. Hence, a 
is the same for different clusters, and it is even almost 
the same for different energy normalizations of the laser 



We have developed a quasiclassical model for small rare 
gas cluster in strong laser fields. This model allows us to 
investigate the influence of several parameters, like the 
atomic element and the cluster size, but also the charac- 
teristica of the applied laser field. We have shown that, as 
a function of pulse length, the energy absorption as well 
as the ionization yield of all but the Ne clusters show a 
clear maximum when the energy content of the pulse is 
kept fixed. This behavior has been attributed to the ex- 
istence of a critical cluster radius R* , whose origin could 
be explained by generalizing the CREI or ENIO concept 
from diatomic molecules to small rare gas clusters. It was 
shown that this mechanism is stable against a change of 
system parameters, even when switching from linear to 
circular polarization. This is a pronounced difference be- 
tween clusters and molecules. 

Finally, we were able to condense the absorption and 
expansion process into a simple model and obtained an 
analytical expression connecting the expansion time and 
the cluster properties. The validity of this expression has 
been confirmed by our numerical data. 

Future investigations will include the transition to the 
plasma regime, which we expect to begin at around 
A" w 10 3 . In this regime, another critical radius should 
emerge, when the laser frequency matches the plasma fre- 
quency. It will be interesting to see how this additional 
critical radius evolves with N and if there is a regime 
where the R* coming from the ENIO process and the 
plasma related R* c coexist. 

We would like the DFG for financial support within 
the Gerhard Hess-program. 
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